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Abstract: Safe and economical disposal of paper mill sludge is a key consideration for forest products industry. A study was conducted to 
examine the effects of amendments of sludge and nutrients on soil surface C0 2 flux (R s ) in northern hardwood forests and to quantify the re¬ 
lationship among R s , soil temperature, and moisture in these stands. The experiment was a randomized complete block design that included 
sludge-amended, fertilized, and control treatments in sugar maple (Acer saccharum Marsh) dominated hardwood forests in the Upper Penin¬ 
sula of Michigan, USA. Results showed that R s was positively correlated to soil temperature (R 2 = 0.80, p < 0.001), but was poorly corre¬ 
lated to soil moisture. Soil moisture positively affected the R s only in the sludge-amended treatment. The Rs was significantly greater in the 
sludge-amended treatment than in the fertilized (p = 0.033) and the control (p = 0.048) treatments. The maximum R s in the sludge-amended 
treatment was 8.8 gmol C0 2 •m' 2 -s' 1 , 91% and 126% greater than those in the fertilized (4.6 pmol C0 2 -m' 2 -s') and control (3.9 |Xmol 
CQ 2 • m' 2 • s' 1 ) treatments, respectively. The Rs did not differ significantly between the fertilized and control treatments. The difference in R s 
between sludge-amended and the other treatments decreased with time following treatment. 
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Introduction 

The pulp and paper industry, a major component of the for¬ 
estry-related economy in the USA, produces large quantities of 
wastewater treatment residuals - or paper mill sludge, which 
must be disposed in an environmentally safe manner. Approxi¬ 
mately 5.3 x 10 6 t of paper mill sludge was produced in 1995 
(NCASI 1999), much of which was disposed in landfills. Many 
of the existing landfills are near capacity, and the cost of creating 
new landfills has increased dramatically due to increased regula¬ 
tions affecting their planning and construction (Camberato et a! 
1997). Therefore, it is a pressing concern for the pulp and paper 
industry to develop alternative methods for disposing paper mill 
sludge that reduce costs without harming the environment. One 
possible alternative to using landfills is the application of paper 
mill sludge to forests. 

Paper mill sludge amendments provide large quantities of 
carbon and nutrients to forest soils and hence affect the decom¬ 
position and mineral cycles. The potential effects of applying 
paper mill sludge on the soil physical and chemical properties 
and leachates have been extensively investigated in forest eco¬ 
systems and agroecosystems (Bockheim et al. 1988; Kraske and 
Fernandez 1993; Simard et al. 1998; Chantigny et al. 1999; Xiao 
et al. 1999; Vance 2000), but the effects of sludge application on 
soil surface C0 2 flux (R s ), the second largest carbon flux in ter- 
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restrial ecosystems (Raich and Schlesinger 1992), are poorly 
understood. 

The objectives of this study were to: (1) examine the effects of 
amendments of sludge and fertilizers on R s in northern hard¬ 
wood forests; and (2) quantify the relationship among R s , soil 
temperature, and moisture in these stands. We hypothesized that 
application of sludge and fertilizers would increase R s because 
increased nutrient availability could stimulate heterotrophic res¬ 
piration (Foster at al. 1980). 

Methods 

Site description and experimental design 

The study was conducted in the Upper Peninsula, near Norway, 
Michigan, USA (45°47'N, 87°45'W). The long-term (29 years) 
mean January and July air temperatures for this area are -11.0 °C 
and 19.8 °C, respectively. The mean annual temperature is 5.2 °C. 
Mean annual rainfall and snowfall are 739 and 1595 mm, respec¬ 
tively. The dominant soil is a mixed typical eutroboralf. 

The forest is mesic northern hardwood, dominated (92% of the 
total density) by sugar maple. Other tree species, in order of 
importance, include American basswood ( Tilia americana L.), 
white ash (Fraxinus americana L.), American beech ( Fagus 
grandifolia Ehrh.), and black cherry (Prunus serotina Ehrh.). The 
understory consists of sugar maple seedlings, as well as assorted 
perennial herbs, with a relatively large presence of leek ( Allium 
tricoccum Ait). General stand structure and soil characteristics 
were summarized in Table 1. 

The experiment was a randomized complete block design with 
three treatments and three blocks. Three blocks were established 
in May 1997, two at the Waucedah (WAUC1 and WAUC2) and 
one at the Nora’s Camp site (NORA2). The three treatments, 
paper mill sludge, fertilization, and control, were randomly as¬ 
signed to one of three plots in each block. The plots at WAUC1 
and WAUC2 were 25m x 25m, and those at NORA2 were 50m x 
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50m. Larger plots were used at the NORA2 site because of larger 
diameter trees and lower tree density (Table 1). 

Table 1. General stand structure and soil characteristics of the study sites 


Characteristics 

WAUCI 

WAUC2 

NORA2 

Geographic location 

45°47'N, 

45°47'N, 

45°47'N, 


87°45'W 

87°45'W 

87°39'W 

Elevation (m) 

315 

315 

320 

Slope (%) 

<5 

15-20 

5-10 

Aspect 

S 

S 

SW 

Stand age (years) 

55 

55 

60 

Stand density (tree - hm ") 

530(65) 

414(155) 

371 (32) 

Basal area (m 2 • hm 2 ) 

20.0(1.5) 

14.9(1.8) 

20.3(1.2) 

PAI (m 2 -mV 

3.6 

2.7 

2.5 

Average DBH (cm) 

25.6(1.1) 

28.4(1.6) 

28.5 (0.5) 

Soil texture 

fine sandy 

fine sandy 

fine sandy 


loam 

loam 

loam 

SOM (%) 

5.1 (0.3) 

7.3 (0.4) 

7.4 (0.7) 

Soil pH 

5.41 (0.13) 

5.74 (0.09) 

5.97(0.18) 

Total N (%) 

0.16(0.02) 

0.08 (0.03) 

0.33 (0.03) 


a. The data inside parentheses are standard errors. WAUC1, WAUC2, and 
NORA2 stand for the research sites at Waucedah site 1,2, and the Nora's 
Camp site, respectively. 

b. PAI = plant area index (one-half plant area index calculated using a 
Li-Cor LAI-2000). 

c. The values for soil are the mean of 9 plots x 4 samples per plot collected 
at 0-20 cm depth. SOM stands for soil organic matter. 

Fertilized and Sludge-amended treatments 

Nutrients were obtained from commercially available fertilizer. 
The “target” fertilizer application rate was to apply N : P : K : S : 
Mg : Ca at a ratio of 100 : 50 : 33 : 9 : 5 : 4 (Gower et al. 1992) 
with N application of 100 kg-hm 2 . However, Ca was applied at 
a slightly higher rate due to the composition of the Mg/Ca fertil¬ 
izer (Gran-U-Lime) (Table 2). 


Table 2. Nutrient application rates in the sludge and fertilization treat¬ 
ments. _ _ 


Nutri¬ 

ent 

Sludge 
(kg • hm' 2 ) 

Fertilizer 
(kg • hm 2 ) 

Fertilizer Source 

c 

3200 

20 


N 

140 

100 

Urea (46-0-0) 

P 

30 

33 

Triple superphosphate (0-46-0) 

K 

5 

50 

(0-0-60) 

S 

50 

9 

Elemental sulfur fertilizer 

Mg 

30 

5 

Gran-U-Lime 

Ca 

770 

6 

Gran-U-Lime 


The application rate for the paper mill sludge was intended to 
provide total N at a rate equal to that of the fertilized treatment 
(100 kg N-hm' 2 ). Samples of paper mill sludge were collected 
from Champion International mill (Norway, Michigan) prior to 
the treatment to calculate the target application rates. However, 
chemical analysis of a paper mill sludge sample obtained imme¬ 
diately after application revealed that the N content of sludge 
amendment was greater than the initial sample (Table 3) and 
39% higher than the fertilization treatment (Table 2). The sludge 
(at a rate of 11.5 t-hrri 2 on dry weight) and fertilizer were ap¬ 
plied on June 17 and June 26 of 1997, respectively. The applica¬ 


tion was done by hand to minimize plot disturbance and ensure a 
homogeneous distribution of sludge. 


Table 3. Nutrient analysis of the paper mill sludge * 


Nutrient 

Concentration 

Nutrient 

Concentration 

C 

28.03 

B 

12 

N 

1.22 

Mn 

869 

P 

0.25 

Fe 

3221 

K 

0.05 

Cu 

10 

Ca 

6.77 

Al 

9914 

Mg 

0.29 

Na 

1935 

S 

0.41 

C/N 

23 

Zn 

1.09 




Note: Nutrient concentrations of C. N, P, K, Ca, Mg, and S are in %, others 
(except C/N ratio) are in mg kg' 1 . 


Soil surface C0 2 flux measurement 

Soil surface C0 2 flux was measured with a Li-Cor 6200 port¬ 
able COs infrared gas analyzer (Li-Cor Inc., Lincoln, NE, USA) 
equipped with a custom 15 cm diameter chamber. R s was meas¬ 
ured approximately every two weeks immediately after the 
treatments, from July 1997 to January 1998 for a total of eleven 
repeated measurement trials that included a large range of soil 
temperatures and encompassed a full growing season. 

Thin-walled (3.2 mm) polyvinyl chloride (PVC) collars (15 
cm inside diameter x 5 cm tall) were beveled on one end to re¬ 
duce compaction and disturbance of the soil and forest floor 
upon placement. Eight collars were inserted 3 cm into the soil at 
random locations in each plot and left in the same locations 
throughout the study for continuity. Rs was measured continu¬ 
ously five times at each collar, with each measurement corre¬ 
sponding to a 5 pmol • m 2 • s' 1 increase in C0 2 concentration 
inside the chamber. Wang et al. (2002) provide a complete de¬ 
scription of the measurement procedures. 

Soil temperature and moisture were measured near the collars 
while R$ was being measured. Soil temperature was measured at 
2 cm (T 2 ) and 10 cm (T 10 ) depths using a digital long stem ther¬ 
mometer (model no. 15-078K. Fisher Scientific, Pittsburgh, PA). 
Soil samples were collected from 0-2 cm and 2-10 cm depths 
near the collars and weighed the fresh weight, and brought back 
to laboratory for soil moisture determination gravimetrically. All 
the soil samples were oven-dried at 70 °C to constant weight and 
weighed to the nearest 0.01 g. 

Statistical analysis 

Statistical analyses were conducted using SAS software ver¬ 
sion 8.0 (SAS Institute Inc. 2000). The mixed effect procedure 
(PROC MIXED) (Littell et al. 1996) was used for the repeated 
measuring analyses, assigning treatment, soil temperature (T 2 or 
T 10 ), soil moisture (W 2 or W ul ), and all possible interactions as 
fixed effects. Block and collar were treated as random effects. 
Natural logarithmic transformation of R s (Ln (R s )) was needed 
to achieve homoscedasticity. A spatial power covariance matrix 
was chosen for the autocorrelation between repeated measure¬ 
ments based on SBC (Schwarz’s Bayesian Criterion). Collars 
within treatment x block were used as the subject in the repeated 
measure analysis. A backward elimination procedure was per¬ 
formed to remove the insignificant (a > 0.05) terms. An analysis 
of covariance (ANCOVA) was conducted to test the overall 
treatment effect on R s by setting temperature and moisture as 
covariates. Added-variable plots (Cook and Weisberg 1999) 
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were used to provide a visual assessment of the effect of soil 
moisture versus temperature on R s . An analysis of variance 
(ANOVA) was conducted to test the treatment effect on R s at a 
specific (measurement trial. 

Results 

Relationships among soil surface C0 2 flux, soil temperature and 
moisture 

T-ratios of the coefficients in the model Ln (R s ) ~ T| () + W m 
were all greater than those in the model Ln (R s ) ~ T 2 + W 2 in 
spite of the similar variability. Therefore, we used T !0 and W| 0 as 
the predictors of R s in this study (The results for T 2 and W 2 was 
not shown). A contrast test (a = 0.05) suggested that the same 
model of R s against T )0 and W 10 was statistically sufficient for 
the fertilized and control treatments, but a model with significant 
different coefficients and terms was required for the 
sludge-amended treatment (Table 4). 


Table 4. Models of soil surface C0 2 flux (Rs) for different treatments *• 

b 


Coefficient 

DF 

Estimate 

SE 

t - ratio 

P 

A 

18 

-0.5719 

0.188 

-3.04 

0.007 

B 

638 

-0.3929 

0.163 

-2.41 

0.016 

C 

15 

0.1071 

0.011 

9.32 

<0.001 

D 

638 

0.0253 

0.006 

4.42 

<0.001 

E 

654 

0.0172 

0.005 

3.80 

<0.001 


a. DF stands for degree of freedom. SE stands for standard error. 


b. The models are: for the control and fertilized treatments: Ln(R s ) = A + 
C*Tio; for the sludge-amended treatment: Ln(R s ) = (A + B) + (C + D)*Tio + 
E*Wio; where Tio and W 10 represent soil temperature and moisture at 10 cm 
depth, respectively. Ln is natural logarithm. 

The Ln (R s ) was positively correlated to soil temperature for 
all the treatments (R 2 = 0.80, n = 84, p < 0.001) (Figure la), but 
was not significantly correlated to soil moisture for the control 
and fertilisation treatments (p = 0.975) (Figure lb). The log (R s ) 
was positively correlated to soil moisture only for the 
sludge-amended treatment (Table 4 and Figure lb). 

Effects of treatments on soil surface C0 2 flux 

The soil surface C0 2 flux for the sludge-amended treatment 
was significantly greater than in the fertilized (p = 0.033) and the 
control (p = 0.048) treatments (Figure 2). The maximum R s for 
the sludge-amended treatment was 8.8 pmol pmol C0 2 • m' 2 • s' 1 , 
91% and 126% greater than those for the fertilized (4.6 pmol 
C0 2 •m' 2 -s' 1 ) and control (3.9 (imol C0 2 •m‘ 2 -s' 1 ) treatments, 
respectively. The R s for the sludge-amended treatment averaged 
2.86 (imol C0 2 • m 2 • s' 1 , and was 38% and 47% greater than for 
the fertilised (2.07 pmol C0 2 -m' 2 -s'') and the control (1.94 
(imol |imol C0 2 ■ m' 2 • s' 1 ) treatments, respectively. However, the 
difference in R s between the sludge-amended treatment and the 
others tended to decrease with the time following the treatments 
(Figure 3). The relative difference in R s between the 
sludge-amended and the control treatments decreased from 116% 
on day 21 to 20% on day 150 following the treatment: and that 
between the sludge-amended and the fertilized treatments de¬ 
creased from 82% to 12%, respectively. There was no significant 
difference in R s between the control and the fertilized treatments 


(p = 0.661). Blocking also had no significant effect on R s (p = 

0.187). 
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Fig. 1 - (a) The natural log of soil surface CO 2 flux (Rs, (imol 
CO • m • 1 • s' 1 ) versus soil temperature at 10 cm depth (Tio, °C). 

The overall statistics are: R" = 0.80. n = 90, p < 0.001. (b) The added-variable 
plot of Ln (Rs) versus soil moisture at 10 cm depth (Wm, %) given the same 
soil temperature at 10 cm depth (Tm, °C). 

Discussion 

Accurate estimates of soil surface C0 2 flux are important for 
constructing stand-level carbon budgets and assessing potential 
impacts of forest management on C0 2 exchange to the atmos¬ 
phere. because R s is the largest C0 2 flux from forest ecosystems 
to the atmosphere (Raich and Schlesinger 1992; Landsberg and 
Gower 1997). This flux is largely the sum of two processes: mi¬ 
crobial respiration associated with decomposition of organic 
matter, and plant root respiration. A number of biophysical and 
biochemical factors affect the two processes differently and in¬ 
teractively. 

Soil surface C0 2 flux is most strongly correlated to soil tem¬ 
perature (Raich and Schlesinger 1992; Lloyd and Taylor 1994; 
Wang et al. 2002). A temperature coefficient (Q !0 ) has com¬ 
monly been used to estimate R s (e.g. Raich and Schlesinger 
1992). Our study also indicated that soil temperature alone ex¬ 
plained 80% of the observed variations in the R s . However, the 
Q 10 model was unable to reflect the significant difference of the 
treatments on R s measured in this study, because the models of 
R s against temperature were statistically identical, i.e. had the 
similar intercepts (p = 0.723) and slopes (p = 0.500), for the three 
treatments (Figure la). Similar results about the inadequacy of a 
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single Qio have been reported in other studies (Howard and 
Howard 1993; Lloyd and Taylor 1994; Rayment and Jarvis 2000; 
Wang et al. 2002). 
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Fig. 2 Average soil surface CO 2 flux (Rs, pmol CO 2 • m" 2 • s' 1 ), soil 
temperature at 2 and 10 cm depth (Tj, T 10 , °C), and corresponding mois¬ 
ture (W 2 , W,„, %) for the sludge-amended, fertilized, and control treat¬ 
ments on the specific sampling day 

The error bars represent the standard errors of the three blocks. Missing stan¬ 
dard errors on Day 12, 1998 is because soil samples were collected for only 1 
block due to heavy snow. ** denotes that Rs for the sludge-amended treat¬ 
ment is significantly (a = 0.05) greater than Rs for both fertilized and control 
treatments; * denotes that R s for the sludge treatment is only significantly (a 
= 0.05) greater than that for the control. 



Days after the Treatments 

Fig. 3 Dynamics of relative differences in soil surface CO 2 flux (Rs, 
pmol CO 2 • m" 2 • s' 1 ) among the treatments following the treatment Bars 

represent block mean (n = 3) with the standard errors. 


Soil moisture is another important environmental control on 
R s , and often interactively influences R s with soil temperature 
(Davidson et al. 1998; Orchard and Cook 1983; Wagai et al. 
1998; Wang et al. 2002). However, this study showed that soil 
moisture did not significantly affect R s in the fertilized and con¬ 
trol treatments unless the moisture was below a critical value 
during the growing season (Figure 2b). When W !0 was less than 
18%, soil moisture exerted a strong limitation on R s for all the 
three treatments. When the W )0 was above this cutoff point, soil 
moisture showed little effect on R s for the fertilized and control 
treatments. For the sludge-amended treatment, R s was positively 
correlated to W, 0 for a broad range (observed from 12 to 38%). 
Although this study was not designed to differentiate the contri¬ 
bution of microbial respiration from autotrophic respiration, we 
did investigate the fine root biomass, the component responsible 
for the majority of autotrophic respiration (Landsberg and Gower 
1997). Fine root net primary production estimated using 
in-growth core technique in 1999, two years after the treatments, 
was statistically similar for the three treatments (Feldkirchner 
2000). We speculated that the Rs difference among the treat¬ 
ments might attribute to the difference in microbial respiration, 
which might have different sensitivity to soil moisture among the 
treatments. 

Soil microbes responsible for decomposing organic matter are 
limited by the availability of carbon substrate (Foster et al. 1980), 
nitrogen (Jackman 1960), and soil acidity (Paul and Clark 1989). 
Paper mill sludge provided a labile carbon source with a low C/N 
ratio (Table 2) that could stimulate heterotrophic activity. We 
observed a substantial increase in R s immediately after the 
sludge-amended treatment, and the difference in R s between the 
sludge-amended and the other treatments declined with the time 
following the treatments (Figure 3). This suggested that the in¬ 
crease in R s flux in the sludge-amended treatment was probably 
attributed to the carbon-enriched sludge, and declined as the 
labile carbon addition was used by heterotrophs, because micro¬ 
organisms first colonize the surface of substrate and process the 
easily available carbon fraction at relatively rapid rate during the 
initial stage of decomposition (Ladd etal. 1996). 

The addition of nitrogen increased soluble carbon and stimu¬ 
lated heterotrophic activity once microbial demand for carbon 
was met (Foster et al. 1980; Jackman 1960). Therefore, it seems 
reasonable that adding nutrients, especially nitrogen, would 
stimulate microbial activity, and hence increase R s . Application 
of ammonium salts often decreased microbial respiration (Fes¬ 
senden et al. 1971; Foster etal. 1980; Roberge 1976; Salonius 
1972). Other studies have reported that R s does not differ be¬ 
tween control and fertilized forest soils (Castro et al. 1994; Fog 
1988; Haynes and Gower 1995; Vose et al. 1995). One possible 
explanation for the similar R s between the control and the fertil¬ 
ized treatments could be inadequate application rate of N, since 
no significant blocking effect in R s was detected even though the 
total N concentration differed by four folds among the three 
blocks prior to the treatments (Table 1). Additionally, soil acidity 
could influence microbial respiration (Paul and Clark 1989). 
However, soil samples taken in the spring of 1998 indicated that 
soil pH was similar among the three treatments (Feldkirchner 
2000 ). 

In summary, soil temperature explained the majority of ob¬ 
served variation in the soil surface C0 2 flux, but was inadequate 
to describe the impact of sludge amendment on C0 2 flux unless 
soil moisture was included in the models. The significant in¬ 
crease in C0 2 exchange to the atmosphere due to sludge amend- 
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ment suggests that caution be taken when considering the appli¬ 
cation of paper mill sludge to the northern hardwood forests 
given concerns on the increase in atmospheric CCL concentra¬ 
tion. 
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